Numerical Simulations of the Invar Effect 
in Fe-Ni, Fe-Pt, and Fe-Pd Ferromagnets 
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The Invar effect in ferromagnetic Fe-Ni, Fe-Pt, and Fe-Pd alloys is investigated theoretically 
by means of a computationally efficient scheme. The procedure can be divided into two stages: 
study of magnetism and calculations of structural properties. In the first stage, an Ising model is 
considered and fractions of Fe moments which point up as a function of temperature are determined. 
In the second stage, density-functional theory calculations are performed to evaluate free energies 
of alloys in partially disordered local moment states as a function of lattice constant for various 
temperatures. Extensive tests of the scheme are carried out by comparing simulation results for 
thermal expansion coefficients of Fei-^Nix with x = 0.35, 0.4, . . . , 0.8, Feo.72Pto.28, and Feo.6sPdo.32 
with measurements. The scheme is found to perform well, at least qualitatively, throughout the 
whole spectrum of test compounds. For example, the significant reduction of the thermal expansion 
coefficient of Fei_ x Ni x as x decreases from 0.55 to 0.35 near room temperature, which was discovered 
by Guillaume, is reliably reproduced. As a result of the overall qualitative agreement between theory 
and experiment, it appears that the Invar effect in Fe-Ni alloys can be investigated within the same 
computational framework as Fe-Pt and Fe-Pd. 

PACS numbers: 65.40.De, 71. 15. Mb, 75.10.Hk, 75.50.Bb 



I. INTRODUCTION 

Fe-based materials are used for various technological 
applications such as springs in watches, car bodies, mag- 
netic cores, and heads of hard disk drives. Despite their 
ubiquity in everyday life, they exhibit intriguing phenom- 
ena that include, among others, high-temperature super- 
conductivity in Fe pnictides 1 , Fermi-liquid breakdown in 
Fe-Nb alloys 2 , and the Invar effect in transition- metal 
alloys 3 . Discovered more than 100 years ago, Invar Fe- 
based materials display anomalously small thermal ex- 
pansion coefficients over broad temperature ranges. Fe- 
Ni alloys with a Ni concentration of about 35 at.% were 
the first to be found 4 . Subsequently, other Invar Fe-based 
materials were reported, some showing ferromagnetism 
(e.g., Feo.68Pdo.32^) and some antiferromagnetism (e.g., 
Fe 2 Ti 6 ). 

Despite the general consensus that the Invar effect in 
Fe-based ferromagnets occurs as a result of magnetism, 
the mechanism giving rise to the Invar phenomenon re- 
mains controversial. Two prominent questions raised by 
recent publications 7,8 have yet to be answered before the 
Invar effect is fully understood: (i) Does the anomaly 
in Fe-Ni appear when changes in the magnitude of local 
magnetic moments with increasing temperature become 
anomalously large? (ii) Are the anomalies observed in 
Fe-Pt, Fe-Pd, and Fe-Ni governed by the same underly- 
ing physics? 

Obviously, any unified theory of thermal expansion in 
Fe-based ferromagnets should capture the Invar effect in 
ferromagnetic disordered face-centered cubic (fee) Fe-Ni, 
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FIG. 1: Average magnetic moment per atom of fee 
Fei-ccNicc at zero temperature plotted against nickel atomic 
concentration, according to the Ising model (circles) and 
experiments^ (triangles). Inset: Concentration dependence 
of calculated Curie temperature (circles) and measured Curie 
temperature^ which has been rescaled by the factor 1.23 (tri- 
angles) . This figure illustrates step 1 of the numerical method 
which we have designed to investigate the Invar effect in fer- 
romagnetic Fe-Ni, Fe-Pt, and Fe-Pd. 



Fe-Pt, and Fe-Pd within a single framework. In princi- 
ple, the linear thermal expansion coefficient of disordered 
fee Fei- X A X with ^4=Ni, Pt, Pd at zero pressure can be 
derived from the Helmholtz free energy which depends 
explicitly on length and temperature. In reality, no ap- 
plications of density-functional theory (DFT) to ab ini- 
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FIG. 2: Estimated fraction of all Fe moments with n Fe first 
neighbors which point up in fee Feo.65Nio.35 plotted against 
n for temperatures below the Curie temperature, Xc. Corre- 
sponding results for Feo.72Pto.28 and Feo.6sPdo.32 are shown 
in insets (a) and (b), respectively. This figure illustrates step 
2 of the method described in Sec. ITTI 
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FIG. 3: Calculated magnetic properties associated with Fe 
sites as a function of temperature in fee Feo.65Nio.35 (black 
lines), Feo.72Pto.28 (blue lines), Feo.6sPdo.32 (red lines), and 
Feo.2Nio.8 (green lines). Fractions of Fe moments which point 
up correspond to solid lines, while demagnetization rates of 
Fe sites are indicated by dashed lines. This figure illustrates 
step 3 of the method. 



tio calculations of finite-temperature free energies have 
been reported to date. One of the major issues in imple- 
menting this strategy is how to incorporate magnetism 
correctly within current approximations to the exchange 
and correlation functional 7,9 . 

In a recent Letter 10 , the magnetic contribution to the 
fractional change in length as a function of temperature 
was studied theoretically for the case of disordered fee Fe- 
Pt. As in our work, the disordered local moment (DLM) 
formalism^ - — was used. However, unlike in our investi- 
gation, effects of lattice vibrations on structural quanti- 
ties were neglected. 

The rest of the paper is organized as follows. First, 
Sec. [Ill introduces a scheme to study the temperature 
dependence of the linear thermal expansion coefficient 
of ferromagnetic disordered fee Fei-^^ with ^4=Ni, 
Pt, Pd. A local- moment model is employed to examine 
magnetic properties; DFT-based calculations and the 
Debye-Griineisen mode l 14 i 15 provide complementary ap- 
proaches for determining contributions to free energies. 
In Sec. IIII| the scheme is tested on alloys with different 
chemical compositions by comparing numerically calcu- 
lated thermal expansion coefficients with experimental 
measurements. Finally, Sec. IIVI summarizes our findings. 
Our work points out the possibility to investigate the 
Invar effect in Fe-Ni, Fe-Pt, and Fe-Pd ferromagnets 
within the same computational framework. 



II. COMPUTATIONAL METHODS 

To study thermal expansion of Fei- X A X over a broad 
temperature interval, we proceed as follows: 

1. We determine the five input parameters which are 
required for the Ising model of the Miiller-Hesse typeife 



the magnitudes of local magnetic moments, Mp e and 
M^, and the three nearest-neighbor exchange constants, 
^FeFe, JpeA, and Jaa- Mp e and Ma are calculated as the 
average magnetic moments on Fe sites and A sites in an 
homogeneous ferromagnetic state by DFT. The method 
employed for J aa depends on the alloying element A. For 
i4=Pt, Pd, Jaa is taken to be zero: this rather crude as- 
sumption embodies the fact that bulk fee metals Pt and 
Pd are both paramagnetic. On the other hand, for ^4=Ni, 
Jaa is taken to be Jaa of ferromagnetic fee N i 17 i 18 . Fur- 
thermore, in all the cases, JFeFe an d J Ye a are determined 
by fitting the calculated zero-temperature average mag- 
netic moment per atom and the calculated Curie temper- 
ature to experimental data^. Before implementing this 
fitting procedure, we rescale the experimental Curie tem- 
perature by a factor of 1.23, to reflect the fact that our 
mean-field solution overestimates the exact Curie tem- 
perature by approximately 23%i£. 

2. Using the above-determined exchange parameters 
and magnetic moment magnitudes, we solve the disor- 
dered Ising model on the fee lattice in the mean-field 
approximation. The disorder is included by using a sep- 
arate mean field for sites with a different nearest-neighbor 
coordination. For example, we allow the mean field to be 
different on Fe sites with 9 Fe and 3 Ni nearest neighbors 
than on Fe sites with 10 Fe and 2 Ni nearest neighbors. 
Thus there are in total 26 mean fields, which are deter- 
mined by a numerical solution of the appropriate self- 
consistency equations. From our solution, we calculate 
the average fraction of Fe moments with n Fe nearest 
neighbors whose local moments are oriented in the 'up' 
direction. 

3. We estimate the fraction of Fe moments which point 
up, XYe-\{T), using the results obtained from step 2. 

4. We perform DFT calculations of the total energy of 
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TABLE I: Theoretical and experimental results for fee Fe-Ni, Fe-Pt, and Fe-Pd. Columns 2 and 3 display fractions of Fe 
moments which point up at several temperatures, according to the Ising model. Columns 4 and 5 show equilibrium lattice 
constants obtained by minimization of total energies (see steps 4 and 5 in Sec.|ll])- Columns 7 and 8 show the result of minimizing 
free energies (see steps 6 and 7 in Sec. [TTJ) . Note that, contrary to column 4, the effect of zero-point lattice vibrations is included 
in the numbers in column 7. Columns 6 and 9 compare the effects of raising XFetOO from its value above the Curie temperature 
to its value at zero temperature on the calculated equilibrium lattice constants ao(xFetOO) an d a(0,XFetOO) f° r eacn an °y 
considered in this table. The rightmost column gives lattice constants measured at 4.2 K 35 i 36 . 
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the random alloy in a collinear magnetic state which re- 
produces the statistics of the local moments' orientations 
xp e ^(T), E (xYe^(T) , a) , for various lattice constants a. 
Depending on the value of XYe\(T), the system is in ho- 
mogeneous ferromagnetic states [case xpetCO = 1]? P ar ~ 
tially disordered local moment (PDLM) states [case 0.5 < 
^FetCO < 1], or DLM states [case x ¥e ^(T) = 0.5]. In the 
two latter cases, up- and down- moments are randomly 
distributed on Fe sites. Total energies are calculated 
using the generalized gradient approximation (GGA)^ 
and within the framework of the exact muffin-tin orbitals 
(EMTO) theory combined with the full charge density 
(FCD) technique—. As in recent theoretical studies on 
Fe-Ni2*2*22 and Fe-P t 1Q i 23 , complete positional disorders 
of chemical species on fee lattice sites and up- and down- 
moments on Fe sites are treated within the coherent po- 
tential approximation (CPA)24. Integration in the irre- 
ducible wedge of the Brillouin zone is carried out over 
several thousands of k-points generated according to the 
Monkhorst-Pack scheme 25 . 

5. We fit the results of step 4 with a Morse function. 
The parameters of the fit give the equilibrium lattice con- 
stant, ao(#Fet(^))> ^ ne bulk modulus, Bq (xpetCO) •> an d 
the Griineisen constant, 7o(xpet(^))- 

6. For each lattice constant chosen in step 4, we 
add to the total energy E{xY^(T),a) a vibrational 



free energy contribution to the Helmholtz free energy, 
F v ib(T, xpetfD? a )- The latter is estimated within the 
Debye- Griineisen model from the outputs of step 5, 
ao(^Fet(T))> Bq{x ¥ ^(T)), and j (x Fe ^(T)). The sum 
of the two terms mentioned above can be written as 

F(T,x Fet (T),a) =E(x Fet (0),a) 
+ F mag (x m (T),a) + F vib (T,x Fet (T),a), (1) 

where 

Fm, g {x F et(T),a) =E(x m (T),a) - £(z Fet (0),a).(2) 

7. We minimize the contribution to the Helmholtz free 
energy (pQ) with respect to a to obtain the equilibrium 
lattice spacing a(T, x Fe ^(T)). 

8. We repeat steps 2 to 7 with different temperatures. 
Subsequently, we apply a cubic-spline interpolation pro- 
cedure. 

9. We evaluate the thermal expansion coefficient 

a(T+ ST, x Fet (T + ST)) - a(T, x m (T)) 
a(T) = hm — -. — -t -, 

v ; st^o a(T,x Fet (T))5T 

(3) 

for the dense set of data from step 8. 
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FIG. 4: Calculated total energy per atom of fee Feo.65Nio.35 in 
a collinear magnetic state E(xp e ^(T), a) plotted against lat- 
tice constant for various temperatures. This figure illustrates 
step 4 of the method. 



FIG. 5: Free energy of fee Feo.65Nio.35 in a collinear magnetic 
state F (T, XFet (T) , a) as determined from Eq. JT]) versus lat- 
tice constant for several temperatures. Black filled circles 
depict equilibrium lattice parameters. 



III. RESULTS AND ANALYSIS 

According to experiments, ferromagnetic Fei-^Ni^ 
with x = 0.35, 0.4, . . . , 0.8, Feo.72Pto.28, and Fe .6sPdo.32 
exhibit a wide variety of thermal behavior, some showing 
the Invar effect ^ 26 ! 27 and others presenting thermal ex- 
pansion similar to that of paramagnetic alloys 26 . For this 
reason, they represent an attractive choice for testing the 
general approach presented in Sec. HU 

We begin the test by considering the input parameters 
of the Ising model. The average magnetic moments on 
each type of site in the homogeneous ferromagnetic bi- 
nary alloys are calculated at zero temperature by means 
of the EMTO method. Our calculated moments on Fe 
sites cover the range from 2.63 /ie for Feo.65Nio.35 to 
2.89 /iB for Feo.68Pdo.32; the moments on Ni, Pt, and 
Pd sites are found to span the interval from 0.3 /hb for 
Feo.68Pdo.32 to 0.64 /iB for Feo.2Nio.8- All these results 
yield fair agreement with available DFT data&i^i 2 ^ 3 -^ 
and experimental measurement s 29 ^ . 

Fig. [1] displays data for Fe-Ni. While the inset com- 
pares calculated Curie temperatures to rescaled experi- 
mental findings, the main panel compares calculated av- 
erage magnetic moments per atom at zero temperature 
to measurements. The quantitative agreement of the 
numerical results with the corresponding experimental 
observations is achieved by varying JFeFe and JFeNi for 
each considered Ni concentration. The fitting procedure 
leads to an enhancement of the Fe-Fe exchange parame- 
ter as the Ni content increases, from a negative value for 
x = 0.35 to a positive value for x = 0.8. This behavior is 
consistent with Monte Carlo simulations 31 . Interestingly, 
when applied to Feo.72Pto.28 and Feo.6sPdo.32, the fitting 
procedure described in step 1 in Sec. [TT] gives ferromag- 
netic coupling between moments on neighboring Fe sites 

(JpeFe > 0). 

Estimated fractions of all Fe moments with n Fe first 
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FIG. 6: Estimated equilibrium lattice parameter of fee 
Feo.65Nio.35 in a collinear magnetic state a(T, XFetf^O) plot- 
ted against temperature T for various values of tempera- 
ture T. The same symbol as in Fig. [5] marks the calcu- 
lated equilibrium lattice parameter a(T, XFetCO) f° r T — 
0, 100, . . . , 1100 K. The thick black solid line results from ap- 
plying a cubic-spline interpolation scheme to a data set which 
is almost twice as large as the number of black filled circles. 
From this curve, we obtain thermal expansion coefficients [see 
Fig. EJa)]. 



neighbors which point up are presented in Fig. [2] for 
n = 0,1,.. .,12; calculated fractions of Fe moments 
which point up are plotted against temperature in Fig. [3j 
Feo.72Pto.28 and Feo.68Pdo.32 are found to exhibit homo- 
geneous ferromagnetism at zero temperature, while the 
magnetic structure of Feo.65Nio.35 appears to consist of 
97% of up-moments and 3% of down-moments. These re- 
sults reproduce available experimental observation s 32 ! 33 . 
For Feo.72Pto.28, Feo.6sPdo.32, and Feo.65Nio.35, the cal- 
culations give XF e ^(T) = 1, 1, 0.9576, respectively, at zero 
temperature. The corresponding values drop by 13, 15, 
and 22% at the reduced temperature T/Tq = 0.75. The 
fraction of Fe moments which point up in the Fe-Ni al- 
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loy is therefore predicted to significantly underestimate 
that of the Fe-Pt and Fe-Pd alloys not only at zero tem- 
perature but also near the Curie temperature. Analy- 
sis of Fig. [2] provides insight into how up- and down- 
moments are distributed among Fe sites for various tem- 
peratures. In the Fe-Ni alloy at zero temperature, down- 
moments are found to reside exclusively on Fe sites with 
11 and 12 Fe nearest neighbors. This picture is consistent 
with recent density-functional total-energy calculations 
performed within the local spin-density approximation 
(LSDA) at the experimental lattice spacing of 3.59 A 8 . 
Perhaps more surprisingly, Fig. [2] reveals that the distri- 
bution of up-moments on Fe sites for any reduced tem- 
perature in the range 0-0.75, more closely resembles a 
random distribution in the Fe-Pt alloy. Accordingly, we 
expect the methodology introduced in Sec. HI] to produce 
more accurate thermal expansion coefficients for Fe-Pt 
than for Fe-Ni. 

In Fig. IH calculated total energy of Feo.65Nio.35 in a 
collinear magnetic state .E^FetCO? °) * s plotted as a 
function of lattice constant for temperature intervals of 
100 K. The estimated values for the equilibrium lattice 
constant ao(#FetOO) a ^ zero temperature and above the 
Curie temperature are reported in Table U along with 
those of other compounds. The curves in Fig. 3] are an- 
alyzed in light of Fig. [3j The equilibrium lattice con- 
stant ao(#FetCO) shifts continuously towards larger val- 
ues with increasing the fraction of Fe moments which 
point up in the system. This confirms expectations 
based on Refs. ^ and |To|- Actually, a behavior simi- 
lar to that observed in Feo.65Nio.35 is seen in each of 
the other systems investigated. To get a rough esti- 
mate of the effect in each alloy, we evaluate the ratio 

Aao/A^Fet = [ao(zFet(0)) ~ a o(0.5)]/[a?Fet(0) - 
The estimated values are displayed in Table HI indicat- 
ing that the dependence on the fraction of Fe moments 
which point up is more pronounced in the Fe-rich alloys 
Feo.65Nio.35, Feo.72Pto.28, and Feo.6sPclo.32 than in the 
Ni-rich alloy Feo.2Nio.8- 

The result of applying the sixth, seventh, and eighth 
steps of the procedure is shown in Figs. [5] and [6] for 
Feo.65Nio.35. In Fig. [5j the total energy £^(xFet(^), a ) * s 
added to the vibrational free energy F v ib(T, x-p e ^(T), a) 
and their sum F(T, xf g ^(T), a) is plotted over a narrow 
range of lattice constants for several temperatures below 
and above the Curie temperature of 614 K. The positions 
of the minima in free-energy curves are marked by black 
filled circles. Results are reported in Fig. [6j In accor- 
dance with experimental observations 34 , it is found that 
the equilibrium lattice constant in the range 0-Tq clearly 
displays a deviation from the monotonically increasing 
behavior seen above the critical temperature. 

While this paper focuses on thermal expansion, it is 
interesting to test whether our simulation properly mod- 
els departures from Vegard's law for Fe-Ni alloys (see, 
e.g., Ref. 35). Similarly to the Invar effect, characteris- 
tic negative deviations from linear behavior in the lat- 
tice constant at very low temperature still await for a 
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FIG. 7: Linear thermal expansion coefficients of fee ferro- 
magnets plotted as a function of temperature. Panel (a): 
Feo.65Nio.35- Panel (b): Feo.72Pto.2s- Panel (c): Feo.6sPdo.32- 
Solid lines show results obtained following the procedure de- 
scribed in Sec. [Ill Dashed lines correspond to experimental 
data 5,27,37 . Vertical arrows indicate Curie temperatures. Our 
observation that these materials all display the Invar effect 
perfectly matches experimental findings. In addition, the the- 
ory correctly predicts the overall trends in a versus T for 
Feo.72Pto.28 and Feo.68Pdo.32. 



complete understanding. As can be seen in Table HI the 
calculated a(0, #Fet(0)) agrees closely with a measured 
lattice constant 35,36 irrespective of the chemical compo- 
sition of the considered material. With the exception 
of Feo.65Nio.35, our data for the Fe-Ni series are fitted 
linearly as a(0,xp e t(0)) versus x. The resulting relative 
deviation from Vegard's law at x = 0.35 amounts to - 
0.14%, which is in quantitative agreement with the tiny 
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FIG. 8: Linear thermal expansion coefficient of fee Fei-xNi^ 
versus nickel atomic concentration, according to the method 
presented in Sec. [Til (filled symbols) and experiments 26 (open 
symbols). Circles and triangles show results for T — 273 K 
and T = 1073 K, respectively. Our approach reliably repro- 
duces the significant reduction of the thermal expansion coef- 
ficient of Fei-ccNix as x decreases from 0.55 to 0.35 near room 
temperature, which was discovered by Guillaume 38 . 



experimental value of -0.08%. 

We now turn to the central question of how well the 
model captures the rich variety of thermal expansion 
phenomena observed in ferromagnetic Fei-^Ni^ with 
x = 0.35, 0.4, . . . , 0.8, Feo.72Pto.28, and Feo.6sPclo.32. 

Figs. [7] and [8] provide a comparison of computed and 
experimentally-determined ^ 27 ! 37 linear thermal expan- 
sion coefficients. While Fig. [7] illustrates the temperature 
dependence of calculated and measured structural prop- 
erties of Feo.65Nio.35 [panel (a)], Feo.72Pto.28 [panel (b)], 
and Feo.68Pclo.32 [panel (c)], Fig. [8] shows how results 
obtained for Fei-^Ni^ vary with nickel concentration at 
fixed temperature. 

The model performs well, at least at a qualitative level, 
throughout the whole spectrum of test compounds. In- 
deed, our observation that Feo.65Nio.35, Feo.72Pto.28, and 
Feo.68Pdo.32 all display anomalously small thermal ex- 
pansion coefficients over broad temperature ranges (i.e., 
the Invar effect) perfectly matches experimental find- 
ings. In addition, the theory correctly predicts the over- 
all trends in a versus T for Feo.72Pto.28 and Feo.6sPdo.32- 
Even the significant reduction of the thermal expansion 
coefficient of Fei-^Ni^ as x decreases from 0.55 to 0.35 
near room temperature, which was discovered by the No- 
bel prize winner Guillaume 38 , is reliably reproduced. 



IV. CONCLUSION 

To investigate theoretically the Invar effect in ferro- 
magnetic disordered fee Fe-^4 with ^4=Ni, Pt, Pd, a 
computationally efficient scheme inspired by previous 
wor k 1Q i 14 has been designed. The procedure can be di- 
vided into two stages: study of magnetism and calcu- 
lations of structural properties. In the first stage, an 
Ising model is considered and fractions of Fe moments 
which point up as a function of temperature are deter- 
mined. In the second stage, DFT calculations are per- 
formed to evaluate free energies of alloys in PDLM and 
DLM states as a function of lattice constant for various 
temperatures. It is worth emphasizing that neither non- 
collinear magnetism 39 nor partial chemical ordering 40 are 
explicitly taken into account at any stage. 

Extensive tests of the approach have been carried 
out by comparing simulation results for thermal expan- 
sion coefficients of Fei-^Ni^ with x = 0.35, 0.4, . . . , 0.8, 
Feo.72Pto.28, and Feo.68Pdo.32 with measurements. De- 
spite a number of approximations (e.g., neglect of static 
ionic displacements^^), the scheme has been found to 
perform well, at least qualitatively, throughout the whole 
spectrum of test compounds. 

As a result of the overall qualitative agreement between 
theory and experiment, it appears that the Invar effect 
in Fe-Ni can be investigated within the same computa- 
tional framework as Fe-Pt and Fe-Pd. This represents 
significant progress compared to previous schemes that 
incorporate DFT calculations. 

In addition, tests results provide evidence that the 
methodology captures the essential physics of the Invar 
effect. For this reason, the present work is currently being 
extended to achieve a better understanding of the phys- 
ical mechanism behind the remarkable phenomenon 42 . 
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